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Abstract. Ultracold quantum gases offer unique possibilities to study interacting 
many-body quantum systems. Probing and manipulating such systems with ever 
increasing degree of control requires novel experimental techniques. Scanning electron 
microscopy is a high resolution technique which can be used for in situ imaging, single 
site addressing in optical lattices and precision density engineering. Here, we review 
recent advances and achievements obtained with this technique and discuss future 
perspectives. 
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1. Introduction 

Scientific breakthroughs has always been associated with technological progress. The 
emergence of more sophisticated measurement techniques enables the discovery of new 
and fascinating phenomena. In particular the exploration of microscopic structures 
greatly benefits from the ever growing capabilities in imaging objects with high spatial 
resolution. At an extremely small length scale, the underlying quantum-mechanical 
nature eventually becomes observable. The realization of Bose-Einstein condensation in 
dilute atomic gases HI El El HI E] and early experiments on it have demonstrated 
that at sufficiently low temperatures the quantum effects are ‘magnified’ such that they 
become observable in real space M. Since then, an explosion of theoretical and 
experimental work has set in, including strongly interacting lattice gases [TOj, interacting 
Fermi gases HD H2] HE], spinor systems HI. artificial gauge fields m HE], quantum 
optical applications im and quantum information processors [181, T9]. 

Obviously, the desire for an in-depth understanding of the numerous phenomena, 
such as superfluidity, quantum correlations, or strongly interacting many-body quantum 
phases drives the quest for novel preparation and imaging techniques. Inspired by 
the principles of scanning electron microscopy, a technique which is unrivaled in its 
spatial resolution compared to optical methods, we have developed an apparatus for 
probing ultracold atomic gases by means of electron-atom scattering [HUSHES]. The 
potential of the technique was demonstrated by the first image of an atomic Bose- 
Einstein condensate with sub-micrometer spatial resolution [22]. Scanning electron 
microscopy of ultracold quantum gases is not only suited for high resolution imaging, 
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but also for the preparation of engineered density distributions in optical lattice sites 
[23] and in vivo studies of many-body quantum dynamics [24, 25]. 

In this review article, we describe how the principles of scanning electron microscopy 
(SEM) can be applied to the manipulation and detection of ultracold atoms and what 
phenomena can be studied with this technique. The structure of the article is as 
follows: In the second section, we give a general motivation for high resolution probing 
of ultracold quantum gases and provide references to the literature for more details. 
Section 3 is of mainly technical nature and introduces the principles of the imaging 
procedure. The relevant electron-atom scattering mechanisms are discussed and the 
working principle of the microscope is explained. We give special emphasis to the 
technical compatibility of SEM when combined with a cold atom experiment. High 
resolution images of ultracold quantum gases and the demonstration of single site 
addressability close this section. Section 4 presents a selection of three experiments 
which highlight different aspects how quantum gases can be studied with SEM: the 
detection of pair correlations, dissipation induced quantum Zeno dynamics and electron 
microscopy of a Rydberg excited Bose-Einstcin condensate. In the last section we discuss 
the perspectives and possible extensions of SEM technique to fermionic quantum gases 
and to the study of open macroscopic systems. 

2. High precision probing and manipulation of ultracold quantum gases 

The realization of Bose-Einstein condensation and Fermi degeneracy in cold gases of 
atoms and molecules has opened a new era in experimental quantum physics [2Bj. 
Experimental and theoretical progress were accompanied by the development of 
detection and manipulation techniques, such as absorption and phase contrast imaging 
or radio frequency spectroscopy [8]. High precision probing and manipulation techniques 
were sought from early on. They are not only indispensable for the study of many-body 
systems on the single atom level, but can also be used to investigate the dynamics, 
the interactions and the coupling to an environment in the quantum regime. In the 
following we give a few examples how experiments can benefit from high performance 
imaging and manipulation techniques. 

One motivation comes from quantum simulation with cold atoms in optical lattices. 
A completely new class of experiments and physical investigations becomes possible 
when atoms in individual sites of an optical lattice can be probed locally. This gives 
access to many characteristic properties of many-body quantum phases [2 T\ such as the 
site occupancy, the tunneling dynamics and the nearest neighbor correlations. Electron 
microscopy [2j3j and high resolution optical microscopy [28] [29] 30] has been developed 
over the past years and has demonstrated the huge potential of this approach. Related 
to these developments there are proposals for quantum information processing , based 
on arrays of single atoms in optical lattices [19]. 

Transport of particles plays a central role in solid state systems. For example, the 
conductivity of a material depends on the quantum-mechanical transport properties of 
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the electrons when moving through the microscopic structure of the material. With 
the help of cold atom systems, fundamental aspects of such transport processes can 
be studied under controlled conditions [an Eg. Tunneling is a common ingredient 
in transport processes and a high resolution imaging technique is ideally suited to 
characterize the microscopic physics and dynamics at a tunneling junction. 

Understanding the non-equilibrium dynamics of closed interacting quantum systems 
lie at the heart of many-body physics [33, 3UIM]- How a closed system transits from an 
initial out-of-equilibrium state to a final equilibrium state is not yet clearly understood. 
The ability to design, control and measure tailored quantum systems on the microscopic 
level paves the way towards understanding local and global properties of non-equilibrium 
dynamics. 

The evolution of an open many-body quantum system is governed by the coupling 
of the system with the environment. During the last decades several advances have 
been made in the study of environmentally induced phenomena like decoherence and 
decoherence-induced selection of preferred states [Ml [37]. Engineering the action of 
an environment as a local dissipative mechanism is an important tool for designing 
many-body quantum state. In situ addressability of quantum gases can therefore help 
to tailor the environmental influence [25], eventually realizing complex decoherence-free 
subspaces and dissipative many-body attractor states. 

The most common methods for high resolution probing and manipulation of 
ultracold quantum gases are based on optical techniques (for details, see [8] and 
references therein). The spatial resolution is diffraction limited and thus depends on the 
numerical aperture of the optical system and the wavelength of the imaging light. For 
current experiments, the highest reported resolution is about 600 nm and is realized by 
quantum gas microscopes [29], (30] • This resolution is sufficient to probe a quantum gas 
with single atom/single lattice site resolution in the Hubbard model regime. 

Many length scales in ultracold quantum gases can even be smaller than the typical 
spacing in an optical lattice, e.g. the healing length, the correlation length, the average 
interparticle distance, the size of a Rydberg wave function or the extent of the Wannier 
function. SEM of ultracold gases provides a technique with a spatial resolution down to 
100 nm which can work in bulk as well as in lattice systems. Moreover, it combines in 
situ and in vivo probing of quantum gases, mass sensitive detection of different atomic 
species and density engineering of ultracold gases. 

3. Scanning electron microscopy of ultracold quantum gases 

This section contains a detailed description of the physical and technical basis of the 
imaging procedure. In the following paragraph, we summarize the most important 
aspects such that the reader can decide to skip the technical part and directly jump to 
the experimental results starting with section 3.3. 

The detection of ultracold atoms by scanning electron microscopy is based on 
electron impact ionization of the trapped atoms with subsequent ion detection (Fig. [I]). 
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Figure 1. Schematics of the working principle, (a) The electron column, represented 
by the dashed line, is based on two magnetic lenses and electrostatic deflection systems 
(see Fig. 3 for more details). It provides a focused electron beam, which intersects the 
atomic sample, prepared in an optical dipole trap, (b) The electron beam is scanned 
across the cloud. Electron impact ioization produces ions, which are guided with an 
ion optical system towards a channeltron detector. The ion signal together with the 
scan pattern is used to compile the image. From Ref. [38] . 


The electron column provides a focused electron beam (EB) with 6 keV electron energy, 
a full-width-at-half-maximum (FWHM) diameter of 100-500 nm and a current of 10- 
500 nA. The smallest beam diameter we could achieve with the column is about 50 nm. 
This size is determined by the electron optics, which are designed to magnify the virtual 
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Figure 2. Time-of-flight spectrum of 87 Rb atoms ionized by electron impact 
ionization. 80% of the produced ions are singly charged. Higher charge states up 
to Rb 7+ are also observed. A peak width of 200 ns ensures that each detected event 
can be assigned to a spatial coordinate of the electron beam. From Ref. jJSj. 


source size of the electron emitter (25 - 30 nm diameter) by a factor of 2.5. This 
magnification was chosen to allow for large beam currents. Different electron optics can 
produce much smaller beam diameters, however at the cost of beam current. The chosen 
configuration is a good compromise between resolution and imaging speed. The electron 
beam is scanned across an ultracold quantum gas of rubidium atoms, which is prepared 
in an optical dipole trap (39]. The ions produced are extracted with an electrostatic field 
and detected by a channeltron [22]. The small diameter of the electron beam ensures a 
high spatial resolution, whereas the ion detection provides single-atom sensitivity. The 
beam is scanned in a TV-like pattern and the temporal ion signal is assigned to the 
beam position in order to construct the image. A typical image sequence takes about 
100 ms, in which a few hundred atoms are detected. The overall detection efficiency 
is limited by the branching ratio between electron impact ionization and non-ionizing 
collisions and the detector efficiency. It amounts to 10 - 20 percent. As the cross-section 
for electron-atom scattering is eight orders of magnitude smaller than the absorption 
cross-section of a resonant photon, the atomic cloud is optically thin for the electron 
beam. For typical parameters, only one out of 500,000 incident electrons undergoes 
a collision. In-depth descriptions of the electron-atom interactions and the technical 
details are provided in the following, typical SEM images are provided in section 3.3. 

3.1. Interaction mechanisms 

In general, the interaction between an electron and a ground state atom can be elastic 
or inelastic. Electron impact ionization of atoms is a special case of inelastic scattering 
which we consider here as an independent process. In an elastic scattering process, 
the scattered atom carries some energy and momentum of the incident electron. The 
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internal state of the target atom remains unchanged but leaves the trap due to the 
transferred momentum. In an inelastic scattering process, the target atom is excited to 
a higher state. The energy transfer to the atom is again accompanied by a momentum 
transfer. In both cases, the atoms escape from the trap and will not be detected as no 
ion is created. In electron impact ionization, the incident electron knocks out one or 
more bound electrons from the target atom leading to the formation of ions, which can 
be detected by an ion detector. It is this process that allows us to detect the atoms. 
Fig. | shows a time-of-flight spectrum of the produced rubidium ions. About 80% of 
the ions are singly charged, whereas the remaining fraction is multiply charged. They 
originate from inner shell ionization events with subsequent autoionization cascades [40]. 

The application of scanning electron microscopy to the detection of ultracold atoms 
has to take into account several constraints: the signal has to be large enough, the time 
resolution must be sufficiently high and the detection process should be efficient in order 
to avoid long measurement times. Multiple scattering events and secondary processes 
might additionally spoil the image and the detected ions have to be assigned to the 
correct position of the electron beam. The initial cold temperature of the atoms and 
the fact that the majority of the ionization processes occurs at small momentum transfer 
to the remaining ion ensure that the created ions have a negligible starting velocity. The 
time-of-flight to the detector is therefore identical for all singly charged ions (see Fig. [2]) 
and the detection time of the ion can be unambiguously assigned to the position of the 
electron beam [38] . 

The signal strength is a crucial parameter. As any scattered atom is removed from 
the gas, the maximum signal recorded at the detector is given by the number of atoms 
in the gas multiplied with the ionization probability. 

The total ionization cross-section for rubidium at 6 keV electron energy is ai on = 
7 ± 0.7 x 10 _1 'cm 2 [44] and represents 40% of all scattering events [421145] . Elastic and 
inelastic electron-atom collisions constitute the remaining events and lead to atom loss 
with no detectable signal. The sum of all three contributions gives the total scattering 
cross section, which amounts to a tot = 1.78±0.14x 10 -16 cm 2 [41]. A detailed description 
of the elastic and inelastic scattering processes and electron impact ionization can be 
found in reference [44] . An atom that is located in the center of the electron beam, has 
a lifetime against collisions with incoming electrons of 


jo&t.ot 

where jo is the current density in the beam center. For typical beam parameters we find 
r = 5 /is, which gives an estimate for a reasonable dwell time per pixel. Most of the 
experiments are performed with a dwell time of 1 or 2 /rs per pixel. This limits the overall 
number of detected atoms but allows for a fast scanning procedure. A scanning speed 
that is faster than the motion of the atoms in a quantum gas (few mm/s) is necessary 
to make sure that the electron beam effectively interacts with an unperturbed cloud. 
Given a full transmission of the ions to the detector and assuming a well adjusted ion 
detector with 50% overall detection efficiency, about 20% of the atoms can be detected. 
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This signal is enough to extract the accurate density distribution from a quantum gas 
but puts limitations on the ability to count atoms. 

In scanning electron microscopy, multiple scattering and secondary scattering 
processes often pose a problem for the proper interpretation of the signal. As the dilute 
atomic cloud collides only with a fraction of 10 -5 of the incident electrons, multiple 
scattering is completely negligible. Concerning secondary processes however, ion-atom 
collisions are of some importance. As the ions are created inside the gas, they can scatter 
with other atoms on their way out of the gas. As the potential between an ion and an 
atom scales as r -4 , where r is the relative distance, the scattering cross section become 
as a large as 10 _12 cm 2 when low temperatures are approached gn sg. For higher 
densities we observe indeed a non-exponential atom loss which indicates the presence 
of secondary collisions [JX]. Depending on the specific imaging mode and investigation 
these processes have to be taken into account. They can be suppressed by reducing the 
transverse extension of the atomic gas or by increasing the electrostatic extraction field 
as the cross section rapidly drops with the ion energy 051 . 

3.2. Experimental setup 

The experiment combines an apparatus for the production of 8 'Rb Bose-Einstein 
condensates with the setup of a scanning electron microscope. We limit the presence of 
magnetic fields that can distort the electron beam by making an all-optical BEC [39] 
and by shielding the main chamber with /i -metal. Except for these additional cares an 
ion detection path is implemented and a test target can be moved into the electron beam 
for calibration purposes. A detailed description of the electron column, the alignment 
and characterization of the electron beam, the ion optics and the detection path is given 
below. 

The electron column. The electron column is mounted on top of the main chamber. 
It is a custom made column with thermal ZrO-Schottky emitter, which provides an 
electron beam with an energy up to 6 keV g6|. The upper part of the column consists 
of three different vacuum chambers: the gun chamber, which holds the electron emitter, 
the aperture chamber, in which a movable stage with several apertures is mounted, 
and the so called intermediate chamber, which consists of a pneumatic isolation valve 
connected to the main vacuum chamber. This valve is open during normal operation, 
and closed when it is necessary to isolate the main chamber from the electron column 
when this one is vented to substitute the emitter (12000 hours lifetime). The lower 
part of the column, consisting of a pole piece made of an iron-nickel alloy extends into 
the main chamber (see Fig. |3|. The current of the electron beam is measured with a 
Faraday cup, which is placed 5 cm underneath the tip. The electron column has two 
magnetic lenses to focus the beam at a working distance of 13 mm below the tip. The 
magnetic field produced inside the second lens, which is close to the atoms, reaches 
2000 G. An efficient shielding provided by the /j -metal pole piece of the tip reduces 
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Figure 3. Technical drawing of the electron column showing various parts required 
for obtaining a highly focused electron beam for the experiment. 


the field to 1 G at the position of the atoms. The effect of spherical aberrations is 
corrected by inserting apertures into the electron beam. Astigmatic aberrations and 
displacements of the electron beam, caused by the lenses, are corrected by electrostatic 
stigmators and deflectors, which are placed in a combined unit behind each lens m- 
The movable aperture strip with apertures of diameter of 300 /irn, 250 /mi, 150 /mi, 100 
/mi, 50 /mi and 20 /irn is placed underneath the first deflector and stigmator stage in 
the aperture chamber together with an electrostatic blanking unit that allows for the 
dumping of the electron beam inside the column with a frequency up to 5 MHz. Two 
physically separated scanning units, consisting of electrostatic quadrupoles, move the 
electron beam over the atomic cloud. The fast scan unit has a field of view of 200 /mi 
x 200 /fin and can be moved with a bandwidth of 10 MHz. The slow scan unit provides 
scans up to 20 kHz and a field of view of 1 mm x 1 mm. 

Alignment and characterization. For the alignment and the characterization of the 
electron beam we use two different test targets mounted on a movable holder. The 
first target is a copper mesh and the second is a hole with a diameter of 200 /irn. A 
conducting plate is placed on the holder to detect secondary and backscattered electrons 
from the test targets. The alignment consists of adjusting the two lenses, deflectors and 
stigmators to maximize the resolution of the image in a convenient field of view. For a 
coarse alignment the mesh target is used and a follow up fine adjustment is performed 
on the hole target. This is also used to determine the beam diameter by taking a set 
of line scans over the edge of the hole. Each of the measured scans is fitted with an 
error function, from which we derive the beam waist. Finally, we remove both targets 
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and optimize the alignment on the atoms in a deep 2D optical lattice with spacing of 
600 nm j46[ ;21, (23j. The highest spatial resolution which we can obtain for a typical 
current of 20 nA is ~ 100 nm. The depth of focus is derived by measuring the FWHM 
of the electron beam at different vertical positions without refocusing. This is achieved 
by moving the hole target along the beam axis by a micrometer positioning translation 
stage. The Rayleigh length is obtained as 35 /mi assuming a Gaussian shape of the 
electron beam with a waist of ~ 100 nm. The vertical extension of the atomic cloud is 
about 6 /irn such that a depth of focus of 35 /mi guaranties a constant electron beam 
diameter across the atomic sample. 

Ion optics and detection. Once the atoms are ionized, they are guided by a series of 
electrodes to a channeltron detector. The ions first hit a conversion electrode (—4.8 kV), 
where secondary electrons are produced and further accelerated towards the channeltron. 
The input of the channeltron has a voltage of -2.2 kV. A detected ion produces a negative 
voltage peak of -10 to -60 mV with a pulse width of about 100 ns. This signal is converted 
into a TTL pulse by a discriminator and amplifier. The pulses are recorded with 10 ns 
time resolution by the control hardware (ADwin-Pro II). 

3.3. Imaging ultracold quantum gases 

A typical image of a Bose-Einstein condensate is shown in Fig. [4j S ummin g over 
many individual images taken in different experimental runs yields a high precision, 
high contrast density profile of the condensate. From the given beam parameters, 
it is obvious that the spatial resolution of a SEM exceeds the diffraction limit of an 
optical microscope. However, an improvement of the resolution is achieved at the cost 
of a smaller beam current. As all characteristic length scales in a quantum gas are 
in the order of a few hundred nanometer (the healing length is ~ 400 nm, the typical 
interatomic distance is ~ 200 nm, the typical spacing in an optical lattice is ~ 500 nm, a 
resolution of 100 nm is sufficient for most purposes. For this reason, we run the electron 
column at working points which are rather unusual for scanning electron microscopes 
and are optimized for the experimental requirements. Typical combinations of beam 
current and beam diameter that are used in the setup are summarized in Table 1. 


beam diameter (FWHM) 

beam current 

90 nm 

12 nA 

250 nm 

100 nA 

400 nm 

180 nA 

5000 nm 

800 nA 


Table 1. Electron beam characteristics for some of the typical settings in the 
experiment. 
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Figure 4. Top: Image of a Bose-Einstein condensate of 10 5 87 Rb atoms. The 
image has 400 x 150 pixels with a pixel size of 300 nm x 300 nm. Each pixel 
was illuminated for 2 gs with the electron beam (140 nm FWHM beam diameter). 
Every dot corresponds to a detected atom. In total, 350 ions were collected during 
the exposure. Bottom: Sum over 300 images. Each image was taken in a separate 
experimental run. From Ref. |22j . 



Figure 5. Image of a Bose Einstein condensate loaded in a two-dimensional optical 
lattice with 600 nm lattice spacing (sum obtained from 260 individual images). Each 
site has a tube-like shape with an extension of 6 /mi perpendicular to the plane of 
projection. The central lattice sites contain about 80 atoms. From Ref. [23] , 


3-4- Single site addressing 

To characterize the resolution of the imaging technique, we load the condensate in a 
two-dimensional optical lattice with 600 nm lattice period (Fig. [5]). The structure of the 
potential is clearly visible and from a quantitative evaluation we can deduce a spatial 
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Figure 6. Engineering arbitrary patterns of atoms loaded in a 2D optical lattice with a 
spacing of 600 nrn. Each emptied site was illuminated with the electron beam between 
1 and 3 milliseconds before the imaging procedure was started. From Ref. [23]. 


resolution of better than 150 nm [39], With such a high resolution, the technique can 
immediately be used for single-site addressing in an optical lattice. To this purpose, 
we selectively remove atoms from individual sites by means of collisions with the 
focused electron beam. In this way, arbitrary patterns of occupied lattice sites can 
be produced [23]j. Four elementary examples are presented in Fig. [6] A single defect in 
the lattice structure is shown in the first panel and represents an ideal starting point to 
study the tunneling dynamics at a single lattice site. The opposite situation corresponds 
to an isolated lattice site and is shown in the second panel. Such a mesoscopic ensemble 
provides, for instance, the possibility to study the transition from few-body systems to 
the thermodynamic limit. It can also act as a paradigm for Rydberg blockade studies 
as the spatial extension of the ensemble is very small. Self trapping and many-body 
tunneling processes can also be studied with such an initial density distribution. A chain 
and a ring of lattice sites are shown in the lower panels in order to illustrate the large 
variety of achievable geometries. The approach allows for any arbitrary pattern that 
fits to the underlying quadratic geometry (Fig. [6|e)). Upon illumination, the number 
of remaining atoms in a site decays exponentially in time. In terms of writing time, a 
lattice site with less than 5 % filling can be prepared in about 1 ms, being sufficient to 
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extensively tailor and structure the quantum gas in the optical lattice. 

4. Experiments 

In this section, we present three examples how quantum gases can be investigated by 
means of SEM. For weak probing (beam current < 50 nA), the electron beam acts as a 
small perturbation to the system and density distributions and pair correlation functions 
can be measured. For strong probing (beam current > 50nA), the back action on the 
atoms becomes important and we observe quantum Zeno dynamics in the system. In the 
third example, we show how the excitation to Rydberg states can even further enhance 
the interaction between the electron beam and a Bose-Einstein condensate. 


4-1. Detection of pair correlations in temporal domain 


The measurement of correlations is an important diagnostic tool to characterize an 
interacting many-body quantum gas. First order correlations are often studied via 
interference experiments. Second order correlations can be determined, e.g., by 
photoassociation [48j, by the detection of intensity noise [J9] or by a single particle 
sensitive detection technique. The second order correlation function g 2 {r ) represents 
the conditional likelihood for detecting a particle at a time r after a previously detected 
particle and quantifies second order coherence. A BEC can be represented as a single 
macroscopic coherent wave function. Hence, for a pure BEC one expects g 2 ( 0) = 1. Well 
above the BEC critical temperature the thermal atoms behave like incoherent particles. 
The bosonic nature of the atoms leads to bunching, g 2 ( 0) = 2. Using scanning electron 
microscopy as a time-resolved local detection method, we can directly measure second 
and even third order temporal correlation functions [21] • To this end, the electron 
beam is pointed continuously in the center of the trapped gas for a given time and the 
time-resolved ion signal is analyzed. 

The general form of the normalized spatio-temporal correlation function of n 
particles at position r, at time U, with i — 1, is given by: 

Mr + + \ (^ t (r 1 ,U)...Tt(r n ,t n )T(r n ,t n )...T(r 1 ,t 1 )) 

g K '(ri,U; ...;r n ,t n ) = --s---—--s--— (2) 

(Tt(n,ti)'F(r 1 ,t 1 ))...(Tt(r n ,t n )T(r n ,t n )) 


where T and HR are the bosonic held operators and (...) indicates the ensemble average. 
We first consider a thermal gas. An analytical expression for (^^(iq, tp, ^,£ 2 ) can be 
derived, extending the approach of Ref. [50] to the temporal domain. For an ideal Bose 
gas at temperature T above the critical temperature T c , trapped in a harmonic potential 
with average trapping frequency lu, we obtain: 


g {1 \r,r) = 


( 1 + A> 


3/2 


exp 


mr 


t c + it 


Thrf 


1 + (U 


(3) 


where r c = j s the correlation time, r = t 2 — t\ and r = r 2 — rj. The above result 
is valid in the limit ujt,Iiuj/( ksT) <C 1. From Eq.[3j we can derive any higher order 






CONTENTS 


14 


1.10 


1 . 08 - 


1.06 


1.04 


1 . 02 - 


1.00 


0.98 



1500 


x (pS) 


Figure 7. Temporal second order correlation function. Data for two different 
temperatures (circles for 45 nK, squares for 100 nK) are plotted together with the 
theoretical prediction. The inset shows the pair correlation function of a BEC. Please 
note that even well below the critical temperature (the thermal fraction cannot be 
detected in time-of-flight absorption imaging) we are able to measure a small residual 
bunching induced by the thermal component. From Ref. |24j . 
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Figure 8. Temporal third order correlation function for a thermal gas. The data 
(dots) shows the result for a thermal cloud with T = 100 nK. The signal has been 
integrated along one time-axis (ti = t 2 — ti, 72 = £3 — ii = 2 ti). The solid line is the 
theoretical prediction. From Ref. [24] . 


correlation function for thermal bosons using Wick’s theorem. In particular, the second 
order correlation can be calculated as g( 2 \r,r) = 1 + \g^\r, t)\ 2 . 

In Fig. [t] we show the measured g( 2 \0,r) := g^ 2 \r) of a thermal cloud of atoms 
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for two different temperatures T = 45 nK and T = 100 nK. The finite extension of the 
cloud along the electron beam axis leads to a reduction of the contrast as the signal is 
integrated over many phase space cells. This can be accounted for by convolving the 
pair correlation function with the density profile of the trapped gas |24j . Fluctuations 
in the total number of detected ions affect the normalization of the correlation function, 
which is obtained by averaging over several experimental cycles. This causes an offset 
shift of the uncorrelated signal to a value 1% above 1. In order to compensate these 
fluctuations we normalize g^ 2 ’ (r) by the factor 1 + cr 2 /(N ) 2 , where cr 2 and (N) are 
respectively the variance and the mean value of the total number of detected ions in the 
different experimental realizations. Experiment and theory are in good agreement and 
the different correlation times are clearly visible. A measurement with a BEC is shown 
in the inset of Fig. [7] and shows no noticable correlations. 

The measurement scheme can also be extended to evaluate third order correlation 
functions. The results are shown in Fig. [ 8 ] The solid line is the theoretical 
prediction based on Eq.[3j g ( 3 ) (ri,T 2 ) = 1 + |g ( 1 ) (ri )| 2 + |g ( 1 ) (r 2 )| 2 + |g ( 1 ) (r 2 - Ti )| 2 + 
2 -Re(g (1) (ti)(7 (1) (t 2 )(7 (1) (t 2 — Ti)), with T] — t 2 — ti and r 2 = t 3 — ti, leaving the global 
amplitude as the only free parameter. As a consequence of the factorial law n\ bunching 
is expected to be more pronounced at higher orders n. For this reason higher order 
correlations can be employed as a highly sensitive test for coherence, with the only 
drawback represented by the need of high statistics. 

Having confirmed the suitability of SEM for the detection of pair correlations, we 
move on to a system with strong quantum correlations: a one-dimensional (ID) Bose 
gas [HI] . We prepare an array of ID Bose gases with help of a blue-detuned two- 
dimensional optical lattice in which we load a BEC. For sufficiently deep lattices, the 
tubes are completely decoupled and the effective interaction strength within each tube 
is given by gm — 2ahuj r , where a is the s-wave scattering length and ay the frequency 
of the radial harmonic confinement. The strength of the correlations within the ID gas 
can be expressed in terms of the Lieb-Liniger parameter 'y(x) = rugm/h 2 p(x), which 
in local density approximation becomes position dependent (x = 0 corresponds to the 
center of the tube). Here, m is the mass of the particle and p(x) is the line density. For 
7 ( 0 ) <C 1 and T below the degeneracy temperature T d = Ti 2 p 2 (0) /2mk B [ 52 ], a weakly 
interacting quasi-condensate phase is formed. For such a phase, the spatial density 
profile in a harmonic trap is expected to be well described by a Thomas-Fermi (TF) 
parabola [53] and g^ 2 \ 0, 0) ~ 1 [52]. In the opposite limit where 7 ( 0 ) 1, the ID gas is 

strongly interacting and can be approximated as a Tonks-Girardeau (TG) gas [M[ [55]. 
For such a gas the density profile is a square root of a parabola and g^ 2 \ 0, 0) <C 1, which 
indicates an effective reduction in the overlap of the particle wavefunctions, resembling 
the fermionic exclusion principle. As a consequence, the local pair correlation function 
is strongly reduced [48], 

We can experimentally access both regimes by varying the density and the radial 
confinement ay. We perform measurements for 7 ( 0 ) ~ 2, representing the onset of 
strong interactions and 7 ( 0 ) ~ 0.5, representing weaker interactions. We measure the 
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Figure 9. Antibunching in a strongly correlated ID Bose gas. Pair correlation function 
for 1 x 10 4 (a) and 6 x 10 4 (b) trapped atoms, representing two different ID interaction 
strengths. The solid lines are an exponential fit to the data. From Ref. Ell- 


pair correlation function in the same way as described above. The results are shown in 
Fig. @ In both cases the antibunching that is due to the repulsive interaction between 
the atoms is visible. The amplitude is only a few percent, as the signal is integrated 
over 8-10 independent tubes. Note that the time-dependent pair correlation function 
for a strongly interacting ID Bose gas cannot be mapped on a corresponding fermionic 
system, because the annihilation operators at two different times do not commute. 

Pair correlations are a sensitive measure for the dynamics of a many-body quantum 
system. We have demonstrated that a local probe such as SEM is well suited to measure 
them. Care has to be taken regarding the prospective measurement time as the low 
detection efficiency enters quadratically in a pair correlation measurement. 

4-2. Dissipation induced quantum Zeno effect 

When the intensity of the electron beam exceeds the weak probing limit, the back action 
on the quantum gas becomes important. Then the electron beam (EB) takes over the 
role of an environment to which the quantum gas is coupled. In the context of the theory 
of open quantum systems, environmental action gives rise to effective Hamiltonians that 
can contain imaginary terms [56l [571 [58] . The losses induced by the electron beam can 
be regarded as a dissipation process leading to such an imaginary potential. As a result, 
a generalized version of the quantum Zeno effect occurs. The quantum Zeno effect can 
be described as the inhibition of the decay of any unstable quantum state by sufficiently 
frequent measurement of that state [59, 60]. If the rate of measurement is high enough, 
the quantum dynamics is slowed down [61] . Here, it is the electron beam that measures 
the position of the atoms in the BEC. 

In the experiment, we locally remove atoms from a Bose-Einstein condensate, 
changing the intensity of the electron beam [25J. The condensate can be described to a 
good approximation by a mean-held wavefunction obeying the Gross-Pitaevskii equation 
(GPE). This remains valid even if the condensate is coupled with the environment. 
Starting from the Lindblad master equation ihd t p = [H,p] + iftCp , where p is the 
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Figure 10. Quantum-Zeno dynamics in a Bose-Einstein condensate. Number of ions 
collected within the first 5 ms of continuous dissipation on a BEC as a function of 
the EB current I. The three panels report the data obtained with w = 127(5), 170(7) 
and 212(8) nm, from left to right. Each data point is the average over 75 experimental 
repetitions. The solid lines are the result from numerical simulations. From Ref. [25). 


density operator of the many-body system, H is the Hamiltonian operator, and £ is the 
dissipation operator such that Cp = — f dx 7 (x)/ 2 ['b + T /5 + pT+T — 2 TpT + ], with y(x) 
the local dissipation rate, we can write the equation of motion for the expectation value 
of the bosonic held operator T as ihd t ( x V) = Tr(A>dtp), which leads to a time-dependent 
GPE with an additional imaginary term 

<90(x,t) 


ih- 


h 2 X7 2 w(x)\ 

+ h ll + 5Wx,i)| 2 -ftVU(x,f). 


(4) 


dt y 2m ' J ' T v " ’ _/l 2 

Here 0 is th e BEC wavefunction, obeying the constraint /|0(x, t)| 2 dx = N(t), V ext 
the trapping potential and g = Anh 2 a/m, a being the s-wave scattering length. The 
technique allows the independent control of the Hamiltonian and the dissipative terms. 

In the experiment, the position of the EB is fixed at the center of the BEC and we 
record the temporal signal from the ion detector. By controlling the beam parameters, 
we can engineer the dissipative term in Eq. [4]as y(x) = la / (2new 2 ) exp(— (x 2 +y 2 )/2w 2 ), 
where I is the EB current, a the electron-atom scattering cross-section, e the elementary 
charge and w the standard deviation of the electron density distribution, assumed to 
be gaussian [25]. In Fig. 10, for three different values of w, we show the number of ions 
collected in the first 5 ms of continuous dissipation as a function of the EB current. 
We observe a non-monotonic dependence of the number of ions produced as a function 
of the EB current for all three dissipation strengths: above a critical value of the EB 
current, the signal decreases despite the stronger dissipation rate. This paradoxical 
behavior is more pronounced for smaller values of w and is the signature of quantum 
Zeno dynamics. Results obtained by numerically solving Eq. [4] additionally taking 
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into account secondary effects like ion-atom collisions [25], are compared in the same 
figure. The very good agreement between the data and the numerically solved model 
demonstrates that the description of the EB as a pure dissipative potential is sufficient 
to capture the main features of the observed dynamics. A localized imaginary potential 
induces total reflection as the strength of the potential U goes to infinity 126. 6 T, !66i 66i. 
The EB current corresponding to the maximum ion production is of special importance, 
since it sets the parameters that make it possible to engineer the most efficient absorbing 
potential. At the same time, the maximum separates two dynamical regimes: for small 
dissipation strength, the unitary Hamiltonian part dominates, while for large dissipation 
strength, the unitary part dominates. 

The results show that the controlled coupling of a quantum gas to an environment 
can lead to interesting dynamics and new ways to engineer many-body quantum states. 
Engineering dissipative processes has already proven to be useful for the generation of 
many-body quantum states ETj. In the future, we expect that such concepts will be 
extended to create quantum states that are not reachable otherwise [68] . 


4-3. Electron microscopy of a Rydberg excited Bose-Einstein condensate 


Rydberg atoms combine strong dipole-dipole interactions with long radiative lifetimes 
and are a promising candidate to implement long-range interactions [69]. The possibility 
to image Rydberg atoms gives direct access to spatial correlation functions and is the 
key component to study new quantum phases m, tailored spin systems nn and energy 
transfer mechanisms [72]. Direct Rydberg atom imaging has been demonstrated so 
far in optical lattices [62], with ion emission microscopy techniques [73], and with a 
scanning autoionization microscope m- Due to the huge electric transition dipole 
moment Rydberg atoms are much more sensitive to electric fields and feature huge 
cross sections for electron scattering. These are therefore well suited for studies in the 
context of SEM. Here, we show that the combination of scanning electron microscopy 
with Rydberg excited atomic samples offers new possibilities to prepare, image, and 
probe theses systems. 

We demonstrate the above mentioned possibilities with a Bose-Einstein condensate 
of rubidium ( N to t = 1.5 x 10 5 atoms) that is off-resonantly excited to the S 8 P 3/2 Rydberg 
state through a one-photon excitation by a high power ultra violet (UV) laser. Details of 
the experimental apparatus can be found in Ref. cst In order to investigate the effect 
of the single-photon Rydberg excitation on the performance of the electron microscopy, 
we perform continuous line scans over the BEC with and without the simultaneous 
illumination by the UV laser. The UV laser is blue-detuned (A = 27 t x 10 MHz) with 
respect to the 5 Si/2 —> S 8 P 3/2 transition. This ensures a long lifetime of the sample and 
avoids the excitation of molecular states consisting of a ground state atom and a Rydberg 
atom that appear on the red-detuned side of the resonance [76] [77]. Illuminating the 
sample with the electron beam and the UV laser, three different single-atom ionization 


channels are possible (see Fig. 11 1 ): (i) direct electron impact ionization of ground state 
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Figure 11. Scanning electron microscopy of a Rydberg excited Bose-Einstein 
condensate, a) Schematics of the three possible ionization channels for a single atom. 
The first process (effective ion rate Ws, green) is direct electron impact ionization of 
ground state atoms. The second ionization channel (effective ion rate Tjjv, blue) is 
a two step process, consisting of an optical excitation to a Rydberg state with Rabi 
frequency Yl followed by a photoionization process with rate 7 opt . The third ionization 
channel (effective ion rate T; c , red) is a three step process that consists of an optical 
excitation to a Rydberg state with Rabi frequency Q, an Z-changing collision (rate "fi c ) 
and electron impact ionization of the Rydberg states (rate 7 eb )- b) Ion rates during 
an electron beam line scan (waist w = 170(20) nm, current / = 25(1) nA) over a BEC 
with N to t = 150k atoms. The red (black) dots show the signal with (without) the 
simultaneous irradiation of the UV laser ( 1 = 7 rnW, w = 700 p, m, A = 2i r x 10 MHz). 
The given rates correspond to the processes described in a) - the tilde indicates that 
they denote the integrated signal coming from all atoms that contribute to the signal. 
From Ref. m- 


atoms, T eb (ii) photo-ionization by the trapping laser following an excitation by the 
UV laser, T uv and (iii) ionization originating from the combined action of the UV laser 
and the EB, ffi c . 

The effect of the individual processes is identified by analyzing the line-scans shown 
in Fig.jllJ). The action of the UV light alone leads to a constant offset in the line scan of 
the Rydberg excited sample and the single-atom ion rate is measured to be Tjjv = 0.3 
Hz. When the electron beam is scanned across the sample and no UV light is present, 
T eb can be directly measured. When both beams are switched on, the additional ion 
rate T/ c appears. As can be seen from fig. 11 :> this leads to an increase of the ion signal 
by a factor of 2. As the direct electron impact ionization of ground state atoms is 
unaltered by the UV light, the increase can only be explained by a much more efficient 
excitation of the Rydberg state by the UV laser. A quantitative analysis m reveals a 
rate of ffi c = 1.4 kHz, which is almost 3 orders of magnitude larger than the excitation 
rate of the UV light alone Yuv = 0.3 Hz. 

In order to understand this giant enhancement of the off-resonant excitation process 
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in the presence of an electron beam we need to look at the interaction between an electron 
and a Rydberg atom. Inelastic collisions of the incoming electrons with the Rydberg 
atom lead to the population of nearby Rydberg states with different angular momentum 
/ and main quantum number n (so-called /-changing collisions). The scattering cross 
section for these processes is huge and is given by [78] 




Here, a 0 is the Bohr radius, m e the electron mass, M the mass of the atoms, Z the 
charge of the nucleus, E the collision energy and R the ionization energy. A(n, /) 
and B(n ,/) are numerical factors that depend on the quantum numbers n and / 
GBj. The collisions can be treated in first Born approximation for high impact 
energies, as present in the experiment. We calculate the /-changing cross section to 
be (T n i{n = 38,/ = 1) = 8.7 x 10 ~ 13 cm 2 for an EB energy of 6 keV. This cross section 
is orders of magnitude larger than the usual electron scattering cross section for ground 
state atoms. The resulting scattering rate into nearby Rydberg states is then given by 


Tic = <rrf/e, 

where e is the electron charge and j 


( 6 ) 

30 A/cm 2 the current density of the electron 


beam. We obtain a scattering rate of 7 j c = 150 MHz, which is much larger than any 
other decay or transition rate of the atom. The dynamics of the excitation process 
is therefore completely dominated by /-changing collisions. This can be cast into an 
effective linewidth given by 7 i c . As the linewidth exceeds the detuning A, the optical 
excitation to the Rydberg state is very efficient. The /-changing cross section is two 
orders of magnitude larger than the electron impact ionization cross section for the 
Rydberg atom, which is given by ai on = 4 x 10 -15 cm 2 mm- Once an atom is excited 
to a Rydberg state, it is scattered several hundred times in neighboring Rydberg states 
before it is eventually ionized. All this happens on a timescale below one microsecond, 
which means that every atom that is excited to a Rydberg state is eventually ionized. 
All three processes together explain the huge enhancement ra- 

The results show that the interaction of electrons with ultracold Rydberg atoms lead 
to drastic effects. The large ionization cross section make such a technique promising 
for the high resolution spatial detection of Rydberg atoms in atomic samples. 

5. Conclusions and future perspectives 

Scanning electron microscopy is a powerful technique for the study of microscopic 
structures and dynamical processes in ultracold quantum gases. It combines superb 
imaging capabilities with novel ways of system preparation. The ability of in situ and 
in vivo probing makes it well suited for the study of local and global non-equilibrium 
dynamics. The finite detection efficiency puts certain practical constraints on atom 
counting applications, but does not constitute a fundamental limitation. Compared 
to state-of-the-art high-resolution optical imaging [291 GO], which features almost 
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hundred percent detection efficiency, scanning electron microscopy provides a number 
of complementary features: (i) the in vivo nature of the technique makes time-resolved 
probing possible. This allows for the measurement of temporal correlation functions - 
an observable that is usually not accessible in globally destructive measurements such as 
optical fluorescence measurements, (ii) The electron beam technique is not restricted to 
two-dimensional arrays of atoms. It is applicable to any species also when fluorescence 
imaging is not possible. Multi-atom multi-orbital physics can be addressed as well, (iii) 
The spatial resolution can in principle be further improved down to a few nanometer. 
The precise shape of on-site wave functions and small-scale spatial correlations can then 
be investigated. 

Scanning electron microscopy can also be applied to fermionic gases, where the 
study of temporal correlation functions, high precision density measurements and non¬ 
equilibrium dynamics help to characterize the different quantum phases. However, the 
magnetic field required for the application of Feshbach resonances can amount to several 
hundred Gauss and the compatibility with the presented electron microscopy technique 
is not obvious. The emerging technology of controlling the s-wave scattering length with 
optical transitions [EDI I8H [82 ] provides an alternative strategy for tuning the interactions 
of a Fermi gas. The recent advent of the dipolar systems erbium [83] and dysprosium 
|83] exhibit low-field Feshbach resonances that are readily compatible with scanning 
electron microscopy. 

A very exciting application of scanning electron microsopy is the study of open 
many-body quantum systems. This includes dissipative many-body models (36] 85] [86], 
m\ and non-Hermitian Hamiltonians [88]. Interestingly, non-Hermitian quantum walks 
have been predicted to show a topological transition [89]. In a more general context, any 
macroscopic quantum system is influenced by the environment and the interplay between 
coherence, decoherence and dissipation determines its properties. For instance, recent 
observations indicate that even biological systems might exploit quantum coherence 
to enhance their functionality [901, 91] . While decoherence is ubiquitous in nature, its 
control is challenging. Scanning electron microscopy of ultracold quantum gases allows 
for a unique way to engineer dissipative processes. This can be used to study generic 
questions of open quantum systems under idealized conditions, such as dark states, 
dissipative attractors and quantum phases in open systems. 
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